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We present a variational Monte Carlo (VMC) method that works equally well for the ground 
and the excited states of a quantum system. The method is based on the minimization of the 
variance of energy, as opposed to the energy itself in standard methods. As a test, it is applied to 
the investigation of the universal spectrum at the van der Waals length scale for two identical Bose 
atoms in a symmetric harmonic trap, with results compared to the basically exact results obtained 
from a multiscale quantum-defect theory. 
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I. INTRODUCTION 

Monte Carlo methods have played an important role in 
our understanding of a variety of quantum systems, espe- 
cially few- and many-body quantum systems with strong 
interactions that are difficult to treat otherwise (see, e.g., 
Refs. [H, 0, H, 13, IE It is also well-known however, 
that most quantum Monte Carlo methods 0, y| are for- 
mulated in such a way that they are strictly applicable 
only to the ground state of a quantum system, a restric- 
tion that has severely limited their applicability. Con- 
sider, for example, the gaseous Bose-Einstcin condensates 
(BEC) of alkali-metal atoms (see, e.g., Q). Any theory 
that intends to treat the real atomic interaction has to 
deal with the fact that the gaseous BEC branch of states 
are in fact highly excited states of a many-atom system. 
There are many branches of states of lower energies, in- 
cluding the first branch of liquid states as suggested and 
studied recently by one of us 

In this paper we present a variational Monte Carlo 
(VMC) method that works the same way for either the 
ground or the excited states of a quantum system. It is 
based on the minimization of the variance of energy, and 
is the method underlying a recent investigation of the 
universal equation of state at the van der Waals length 
scale [1, |]| for few atoms in a trap ■ The details of 
the method were skipped in the earlier article [10l |. both 
because the focus there was on a single gaseous BEC 
state, which was not the best example illustrating the 
method, and because there were no other independent 
results to directly compared with, except in the shape- 
independent limit [10 ] . 

We present here, in Sec. II, the details of the vari- 
ational Monte Carlo method based on the minimization 
of the variance of energy and shows that it applies equally 
well to the ground and the excited states of a quantum 
system. In Sec. Ill, we present a better illustration of 
the method through the universal spectrum at the van 
der Waals length scale for two identical Bose atoms in a 
symmetric harmonic trap. It is an example where results 



for multiple energy levels can be obtained independently 
using other methods [ll|, [l^, [H, EH , including, in partic- 
ular, a multiscale quantum-defect theory (QDT) p^lflq. 
Conclusions are given in Sec. |lVl We point out that in 
the process of writing this article, we have discovered 
that an equivalent approach has been developed earlier 
by Umrigar et al. 17]. The derivation of our method, 
and the applications presented here and earlier [lOj, are 
however different. 



II. VARIANCE MINIMIZATION VARIATIONAL 
MONTE CARLO METHOD 

Consider the time-independent Schrodinger equation 
H\9 n )= E n |¥„) , (1) 

where the energy eigenstates \^ n ) form a complete, or- 
thonormal basis. 

Existing quantum Monte Carlo methods are mostly 
based on the fact that for an arbitrary trial wave function 
satisfying proper boundary conditions, we have 
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which means that the ground state wave function is the 
one that minimizes the energy functional Et[9t]- The 
proof can be found in standard quantum mechanics text- 
books (see, e.g., [Hj]). 

The variance minimization variational Monte Carlo 
method (VMVMC), as proposed here, is based on the 
functional 
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The proof of Eq. ([3]) and its physical meaning can be best 
understood by expanding the trial wave function using 
the complete basis defined by Eq. |T]) to write 77 \^t\ as 
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From Eq. ([4]), it is clear that zero is the minimum of 
the functional jj^t]; an d this minimum is reached when 
and only when Et — E n and (^ m \^fT) = for m ^ n, 
namely, only when \^St) is an eigenstate of energy as de- 
fined by Eq. ([1]). This statement is equally applicable to 
the ground and the excited states of a quantum system. 

The implementation of VMVMC, based on the mini- 
mization of the variance of energy ^['J't], is straightfor- 
ward. It does not require much more than the standard 
VMC, as we illustrate here using the example of identical 
particles. 

Consider N identical particles in an external potential 
and interacting via pairwise interactions. It is described 
by a Hamiltonian: 



N 



N 



with 



i<j=l 



hi = ~ — Vf + Kartfo) 



(5) 



(6) 



Here V ex t(r) is the external "trapping" potential, and 
v(r) is the interaction between particles. 

For the evaluation of the energy functional, we have 



(^t\H\^ t ) = (^T\Nhi + -N(N - ]>i 2 |*t) 

dT*J* T — {Nhi + -N(N - l)w(n 2 )}*T 



JdT^ T E Loc {T) , 



(7) 



where r represents an N particle configuration specified 
by their 3N coordinates. Eloc is the so-called local en- 
ergy, and is given by 



} . ^V?* T +^VV ext (r 1 )+iiV(iV-l) U (ri2) 
2m / Wt 2 



The average energy is therefore 

/ dT^ T E Loc (T) 



En 



(8) 
(9) 



This is the standard integral in VMC, and can be eval- 
uated using standard Monte Carlo methods such the 
Metropolis method (see, e.g., [Hj]). 

In order to calculate the variance of energy, one must 
also determine the average of H 2 . This can be done by 
first noting that, similar to Eq. (|T|). we have 



* m |ff|* T ) = / dT^ T E Loc {r) , (10) 



where |$m) is an eigenstate of energy as defined by 
Eq. (fT]). We have therefore 



y m \H\y T ) (y m \H\y T 



= E / drdr'[^ m (r'W T (r')E* Loc (r') 

m 

X ** m (T)* T (T)E Loc ( T )} . (11) 

Using the completeness relation 

J2^m(r'W m (r) = S(r'-T), (12) 



dT^* T ( T )^ T {T)\E Loc (T)\ 2 , (13) 



we obtain 

and therefore 

'* T \H 2 \* T ) j dTn{T) y T{T)lELoc{T)l 2 
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The computation of the variance of energy, Eq. (J3]) , has 
thus been reduced to two integrals, Eqs. (J7]) and (fl4|). 
both of which involving the same local energy, El oc , that 
one encounters in standard VMC. It is clear that the for- 
mulation and the equations in this section are applicable 
to both bosons and fermions. 

One can easily show that our method is equivalent to 
that of Umrigar et al. [17]. However, we believe that 
our derivation provides a more rigorous foundation and 
shows more explicitly why it works for both the ground 
and the excited states. 



III. SAMPLE RESULTS FOR IDENTICAL BOSE 
ATOMS IN A SYMMETRIC HARMONIC TRAP 

The VMVMC, as outlined in Sec.HH was first applied 
in Ref. [l0| to study the universal equation of state at 
the van der Waals length scale [8|, [9( for few identical 
Bose atoms (N — 3-5) in a trap. To better illustrate 
and to further test the method, we investigate here the 
universal spectrum at the van der Waals length scale for 
two identical Bose atoms in a symmetric harmonic trap. 
It is a problem for which accurate results can be obtained 
independently using a variety of methods [TJ E EH EH j 
including a multiscalc QDT [l5|, [l(| . 

Two identical Bose atoms in a symmetric harmonic 
trap are described by the Hamiltonian, Eqs. ©-flS]), with 
N = 2, and 



Vextiji) = -mij 2 r 2 , 



(15) 



where m is the mass of an atom, and uj is the trap fre- 
quency. 

For the trap states of interest here, we take the trial 
wave function to be of the form of 



[0i(ri)0 2 (r2) + 0i(r 2 )^ 2 (ri)] F{r 12 ) 



(16) 
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where <j>i and <f>2 are independent-particle orbitals, and 
F is the atom-atom correlation function that is discussed 
in more detail in Ref. [13]. Specifically, we use 



F(r) 



Au\(r)/r , r < d 
r> d 



(r/d) 7 

where u(r) satisfies the Schrodinger equation: 



m dr 



A 



u\(r) = , 



(17) 



(18) 



for r < d. 7 is the parameter characterizing the long- 
range correlation between atoms in a trap, with 7 = 
(meaning F = 1 for r > d) corresponding to no long- 
range correlation. Both d and 7 are taken to be varia- 
tional parameters, in addition to the variational param- 
eters associated with the descriptions of 4>\ and 02 . The 
parameters A and A are not independent. They are deter- 
mined by matching F and its derivative at d. Our choice 
of F differs from traditional choices (see, e.g. Ref. @) not 
only in its treatment of the short-range correlation, but 
especially in its allowance for the long-range correlation 
characterized by parameter 7. This was first suggested 
by a multiscale QDT treatment of two atoms in a sym- 
metric harmonic trap [lol [l6j , and was later found to be 
the key for treating N trapped atoms in cases of strong 
coupling, namely when the s wave scattering length clq 
becomes comparable to or greater than the trap length 
scale a,ho — (fi/tnuj) 1 / 2 [Toj . 

For atoms in their ground state, the atom-atom inter- 
action is of the van der Waals type of —C n /r n with n = 6 
at large interatomic separations, i.e., 



v(r) 



-Cth 



(19) 



This interaction has an associated length scale of (3§ = 
(mCe/h 2 ) 1 ^, and a corresponding energy scale of se = 
(h 2 /m)(l//3 6 ) 2 [13]. Over a wide range of energies that is 
hundreds of se around the threshold [2l|, [22j , the details 
of atomic interactions of shorter range than /3g are not 
important, and can be characterized by a single parame- 
ter that can be the s wave scattering length ao, the short 
ran ge K matrix K c , or some other related parameters 
[2H |23|. H3]. In this range of energies, the spectrum of 
two atoms in a trap follows a universal property that can 
be characterized by [T3, [13, E3] 
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fli(a /a ho , f3 6 /a h o) , 



(20) 



and is called the universal spectrum at length scale (3§. 
Here VLi are universal functions that are uniquely deter- 
mined by the number of particles, the exponent of the 
van der Waals interaction (n = 6), and the exponent 
of the trapping potential (2 for the harmonic trap). The 
strengths of interactions, characterized by Cq and to, play 
a role only through scaling parameters such as /3@ and 

dho- 
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FIG. 1: The universal spectrum at length scale f3e for two 
Bose atoms in a symmetric harmonic trap as a function of 
ao/iho for 06 / dho = 0.001. Solid line: results from a multi- 
scale QDT [3|. Symbols: results of VMVMC. 



As in Ref. [IJJ], the universal spectrum at length scale 
/?6, namely the fij's in Eq. (j2TJ)) . can be computed by using 
a correlation function, Eq. (fTF]) , with u\(r) as given by 
the angular-momentum-insensitive quantum-defect the- 
ory (AQDT) 0, 



u Xs (r s ) = B[f x %(r s ) - K c gZL (r s )} 



c(6) 



(21) 



Here B is a normalization constant, ffflp and j are 
universal AQDT reference functions for —Cs/r e type of 
potentials [3, H3]- They depend on r only through a 
scaled radius r s = r//3 6 , and on energy only through a 
scaled energy X s = X/se- K c is the short-range K matrix 
21 1 that is related to the s wave scattering length a by 

24, m 



j*r(i-&) 



r(i + 6) 



K c + tan(7r&/2) 
K c - tan(?r&/2) 



(22) 



where b — l/(n — 2), with n = 6. 

Figure [T] shows a portion of the universal spectrum at 
length scale /3q for two Bose atoms in a symmetric har- 
monic trap. Specifically, it gives the energies of the first 
three s wave trap states as a function of ao /dho- The cor- 
responding fas used in Eq. (|16p are independent-particle 
orbitals based on standard solutions for a single particle 
in a symmetric harmonic potential (see, e.g., [261]). For 
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FIG. 2: The same as Fig. [Qexcept for f3§/a ho = 0.1. 

the lowest s wave trap state, they are taken to be 

&(r) =exp{-a l x 2 ) , i= 1,2. (23) 
They are taken to be 

4>i(r) = exp(—aix 2 ) , 

02 (r) = Q - x 2 ^) ea;p(-a 2 x 2 ) , (24) 

for the first excited s wave trap state, and 

M T ) = (| - ea^-OjSC 2 ) , i = 1,2, (25) 

for the second excited s wave trap state. Here a; is a 
scaled radius defined by x — r/aho- The variational pa- 
rameters are d, 7, a\, and «2 in all three cases. The 
variance of energy is calculated according to Sec. HH and 
the minimization is carried out using a type of genetic 
algorithm. 

Both Figs. □ and H show that the results of VMVMC 
are in excellent agreements with those of a multiscale 

qdt mm. 

which gives basically exact results for two 
atoms in a symmetric harmonic trap. (The scaled energy 
per particle, £j/(27ki>), used here is related to the scaled 
center-of-mass energy, e = e/hui, used in Ref. [ll|, by 
Ei/(2huj) — (e + 3/2)/2.) The agreements are all within 
the variances of energy, which are smaller for weaker cou- 
pling (smaller ao/aho) and greater for stronger coupling, 
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FIG. 3: A comparison of the spectra for two different values 
of l3§/aho, illustrating the shape-dependent correction that 
becomes more important for greater values of f3§/aho and for 
more highly excited states. 



but are in any case less than 1.8 x 10~ 3 for all parameters 
considered. The results shown in Figure [TJ which are for 
a small fle/aho = 0.001, illustrate the shape-independent 
limit of Pe/dho — > for states with Ei/2 ~ Ticu <C Se 
[ToL Il6j |. They agree, in this limit, with the results ob- 
tained using a delta- function pseudopotential 11]. For 
greater Pe/dho, the effects of the van der Waals inter- 
action become gradually more important, especially for 
strong coupling (ao/aho ~ 1 or greater) and for more 
highly excited states [13, EH . This is illustrated in Fig- 
ure [31 which compares the results for Pe/ai lo = 0.1 with 
those for j3§/ah = 0.001. We note that even the low- 
est trap state is itself a highly excited diatomic state. 
There are other "molecular" states that are lower in en- 
ergy [lj|[l6j]. This fact does not, however, lead to any dif- 
ficulties because VMVMC works the same for the ground 
and the excited states. It is for the same reason that we 
were able to investigate the gaseous BEC state for few 
atoms in a trap (Io| . which is again a highly excited state. 
More detailed discussions of the universal spectrum at 
length scale (3q for two atoms in a symmetric harmonic 
trap, including the molecular states and the spectra for 
nonzero partial waves, can be found elsewhere [161 ]. 
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IV. CONCLUSIONS 

We have presented a variational Monte Carlo method, 
VMVMC, that works the same for the excited states as 
it does for the ground state. The method is tested here 
through the universal spectrum at length scale (3q for two 
identical Bose atoms in a symmetry harmonic trap, for 
which the results from VMVMC are found to be in excel- 



lent agreements with the basically exact results derived 
independently from a multiscale QDT [TH, [l6| . 
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